Snap-Through of Unsymmetric Laminates Using Piezocomposite Actuators 
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ABSTRACT 


The paper discusses the concept of using a piezoceramic actuator bonded to one side of a two- 
layer unsymmetric cross-ply [0/90] T laminate to provide the moments necessary to snap the laminate 
from one stable equilibrium shape to another. This concept could be applied to the moiphing of 
structures. A model of this concept, which is based on the Rayleigh-Ritz technique and the use of 
energy and variational methods, is developed. The experimental phase of the study is discussed, 
including the measurement of the voltage level needed to snap the laminate. The voltage 
measurements and shapes are compared with predictions of the models and the agreement between 
measurements and the predictions are reasonable, both qualitatively and quantitatively. Suggestions 
for future activities are presented. 
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INTRODUCTION 

As is well known, a thin, unsymmetrically laminated, elevated temperature cure, fiber-reinforced 
composite laminate with no external loads, such as a two-layer cross-ply graphite-epoxy [0/90] T 
laminate, can have multiple equilibrium shapes when cooled from the curing temperature to a lower 
operating temperature; that is, the laminate can be multistable. In the case of a cross-ply laminate, 
when the geometry is such that the laminate is multistable, there are two stable cylindrical shapes and 
one unstable saddle shape. The laminate can be changed from one stable cylindrical shape to the other 
by a simple snap-through action by applying moments along opposite edges of the laminate. The 
unstable shape, of course, is never achievable. This multistable phenomenon has been studied by a 
number of investigators [1-4], and is a result of the coupling of the residual stresses due to cooling 
with the geometric nonlinearities due to the large out-of-plane deflections involved. 

This study considers the use of the NASA Langley Research Center Macro-Fiber Composite 
Actuator (MFC), an orthotropic piezoceramic actuator, to effect snap-through behavior of a [0/90] T 
laminate; MFC actuators are prepackaged and the actuation is accomplished by applying a voltage to 
the actuator. The multistable characteristic could be applied to the design of morphing structures. The 
geometry of the structure that was considered in this study is shown schematically in fig. 1; the MFC 
actuator is centrally bonded to one side of the laminate to form the laminate-actuator combination. 
Predictions from a model of the laminate-actuator combination based on the Rayleigh-Ritz technique 
will be compared with experiments. In particular, predictions of the shapes of the laminate-actuator 
combination and of the voltage applied to the actuator that causes snap-through are compared with 
experiments. 
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THE MACRO-FIBER COMPOSITE ACTUATOR 


The MFC actuator consists of piezoceramic PZT macrofibers embedded in a structural epoxy 
matrix and sandwiched between interdigitally-electroded polyimide films [5]. Two-dimensional 
modeling of the piezoelectrically-induced dilatational strains of the MFC actuators used in this study 
can be accomplished in a manner similar to the typical modeling of a monolithic piezoelectric wafer, 
namely, 
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where dn and d u are the effective piezoelectric constants for the MFC actuator, AV is the voltage 
applied to the MFC actuator, and Ax, is the electrode spacing. The term effective is used to describe 
the piezoelectric constants because the interdigitated electrodes create a nonuniform electric field in 
the actuator, but by using AV/Ax t , an average, or effective, electric field in the actuator is considered; 
as with properties of other composites, the electomechanical properties of the MFC actuator are 
“smeared” within the volume of the actuator. In the above strain notation, and throughout this paper, 
the principal material coordinate system nomenclature standard to fiber-reinforced composite 
materials is used. That is, the material coordinate system is denoted by the subscripts 1, 2, and 3, and 
the global coordinate system is denoted by x, y, and z, as shown in fig. 1 . 


MODEL DEVELOPMENT 

The analytical model that was developed for this study was divided into three parts: cooling the 
cured laminate, bonding the actuator to the laminate, and applying voltage to the laminate- actuator 
combination. As in [2], in accordance with the requirements of the Rayleigh-Ritz technique, the 
current model relied on using good approximations of the displacement fields in conjunction with an 
extension of classical lamination theory by including geometric nonlinearities in the strain- 
displacement relations. Stationary values of the total potential energy of the laminate are sought. It 
was believed that a ten-parameter approach could be used for predicting the shapes of the laminate- 
actuator combination. Accordingly, in each of the three steps, the three components of displacement 
of the geometric midsurface of the laminate, the reference surface in this analysis, were approximated 
by 
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( 2 ) 



where u°, v°, and w° are the displacements in the x-, y-, and z-dircctions, respectively. The ten 
parameters c\ through c 10 are unknown, but to-be-determined, coefficients. Because of symmetry 
requirements on the displacement, the third-order vv' -displacement polynomial could not be used to 
represent the 11 / -displacement field over the entire laminate. However, a quarter- symmetry argument 
was used to allow the use of this polynomial to represent the vc' -displacement field, and the total 
potential energy was calculated over only one-quarter of the laminate. It should also be mentioned 
that, consistent with eq. 2, the laminate was considered to be fixed at the geometric center, and with 
no external applied loads. 


Step I: Laminate cooling from cure 

The first step of the model was to determine the initial cooled shape of the two-layer [0/90] T 
laminate. Because this step would not include the MFC actuator, and no forces were applied to the 
laminate, only thermal effects within the laminate needed to be considered. Accordingly, the total 
potential energy, II, used in Step 1 was the strain energy of the cooled laminate, IT!, and is given by 

(l A A 3 

n = n 1 = 4 2 J 0 2 it, 2 i- 2 [ (Cr ^ _ a V K + iPy - o]. ) £ y + iPxy “ G ly )Y xy ]dxdvdz (3) 

V 

where the as are the stresses in the laminate, the a T s are the so-called thermally-induced stresses, the 
ss and y xy are the strains, and L x and L v are the sidelengths of the laminate when flat as depicted in fig. 
1 . Also as shown in fig. 1 , the thickness coordinates z 0 and z 2 are the coordinates for the bottom of the 
0° graphite-epoxy layer and the top of the 90° graphite-epoxy layer, respectively. 

The stresses in eq. 3 are given by 

= 011^+012^+016^ ~°l 
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where the Q s are the transformed reduced stiffnesses of the graphite-epoxy layers and the thermally- 
induced stresses are given by 
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r xy - oc xy AT = 0 

The as are the coefficients of thermal deformation of the layer and AT is the temperature change from 
the cure temperature to room temperature, here the operating temperature. Note that Q l6 , Q 2( , and a xv 
are retained in the formulation, even though they are zero for 0 and 90° layers. 

For small strains and moderate rotations, the midsurface strains and curvatures in the laminate are 
given by the von Karman approximation to the more general strain-displacement relations, namely, 
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According to the Kirchhoff hypothesis, the strains as a function of the thickness location, z, are given 
by 
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By using eq. 2 in eqs. 7 and 8 and, in turn, substituting these results into eqs. 4 and 3, taking into 
account eqs. 5 and 6, and carrying out the integrations of eq. 3, the total potential energy is reduced to 
an algebraic equation in terms of material properties, geometry, and the coefficients c\ through c\o. 
These coefficients are determined by solving the ten simultaneous nonlinear algebraic equations that 
result from equating to zero the first variation of total potential energy with respect to these 
coefficients, namely, 
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Stability of the solution is determined by checking the positive definiteness of the ten-by-ten matrix 
associated with the second variation of the total potential energy with respect to these ten coefficients. 
The algebraic manipulations, integrations, and differentiations in this and the following steps were all 
accomplished with the aid of Mathematica [6]. 


Step II: Bonding the MFC actuator to the laminate 

The MFC actuator was to be bonded to the laminate by placing the laminate and actuator, with 
adhesive, inside a vacuum bag, and using a vacuum pump to evacuate the bag to pull the laminate and 
the actuator into contact until the adhesive cured. The result was that the contacting surfaces of the 
actuator and the laminate would develop the same curvature during the bonding process, the laminate 
losing some curvature and the actuator gaining curvature. This common curvature would initially be 
unknown. Step II of the model was developed to simulate this vacuum bonding technique, i.e., to 
determine the shape of the laminate-actuator combination resulting from the vacuum bonding process. 
This resulting shape of the laminate-actuator combination when the adhesive had cured will be 
referred to as the first actuator-added shape. With the actuator bonded to the laminate, there would be 
another unknown shape to which the laminate- actuator combination could be snapped. This other 
shape will be found in Step III, and will be referred to as the second actuator-added shape. 

The total potential energy of the laminate and MFC actuator bonded to it in Step II was the sum 
of the strain energy from eq. 3, Tf, and the strain energy of the actuator, which is given as 


n 2= 4 ^{ 0 2 | 0 2 ^\o?el 2 +afef +a a *y%yxdydz 


where cr “ 2 , a a v 2 , and rr “ 2 and s a x 2 , sf , and y a2 are the stresses and strains in the MFC actuator 
due to the bonding process. As shown in fig. 1, the inplane dimensions of the actuator when flat are 
L w< and L\" ( . The thickness coordinates z 2 and z 3 define the thickness of the actuator. Until the 

adhesive began to cure, the actuator could slip relative to the laminate, so the stresses created in the 
actuator due to the bonding process would be just bending stresses. The stress-strain relations for the 
bonded actuator are 
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For simplification, it is assumed that the curvatures of the midsurface of the laminate and the 
midsurface of the actuator (rather than the curvatures of the contacting surfaces) are the same. The 
strains in the actuator may then be approximated as 



where the ks are not known, but will be the same in the actuator and the laminate. The displacement 
fields of the form of eq. 2 are again assumed to be valid. For Step 11, the total potential energy of the 
laminate-actuator combination is given as 

n = n, +n 2 (13) 

It should be noted in eq. 1 3 that, in the total potential energies of the laminate and actuator, the work 
terms due to the stresses between the laminate and actuator cancel each other. The coefficients c\ 
through c io for this step are again determined by finding stationary values of n, as in eq. 9, and 
stability is checked by examining the second variation. 

It is important to again delineate at this point the consequences of bonding the actuator to the 
laminate. The analysis shows that when cooled, the two-layer [0/90] T cross-ply laminate considered 
could have either one stable equilibrium shape, or it could have an unstable saddle shape and two 
stable cylindrical shapes, depending on the dimensions of the laminate. In the latter case, the two 
cylindrical shapes would have equal and opposite radii of curvature and perpendicular generating 
axes. The actuator could be bonded to the cooled laminate when the laminate was in either of these 
two cylindrical shapes. As a result of bonding the actuator to the laminate with the laminate in a 
particular shape, that particular shape would change, and the shape to which it could be snapped to 
would also change. The shape resulting from Step II, the first actuator-added shape, would be the 
initial shape for Step III. 

Step III: Actuation of the MFC 

Step III of the model considered the laminate- actuator combination with no slipping between the 
laminate and the actuator, and included the deformation due to voltage applied to the actuator. The 
initial shape was to be specified by the coefficients c\ through cio associated with the shape 
determined from Step II; for Step III, these known coefficients were renamed c [ through cj 0 , 
respectively. The superscript ‘f denotes ‘initial.’ The total potential energy of the model was again 
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broken into the laminate contribution, again Hi, and the actuator contribution, n 3 . The contribution to 
the total potential energy from the MFC actuator in Step 111 is written as 


n 3 = 4| 
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where the s ai s and y " 3 are the strains in the actuator and are to be discussed below. The midsurface 

of the graphite-epoxy laminate remains the reference surface. The <r fl3 s are the stresses in the 
actuator and are given by 
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The a E s are piezoelectrically-induced stresses and are given by 
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where the ds and y 2 are piezoelectrically-induced strains. The strains in the actuator, <s'" 3 s and y a xy . 
are given by 


„a 3 
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where the ss and y xy are Kirchhoff-like strains that are zero at the reference surface (i.e., at the 
midplane of the laminate) and vary linearly with thickness. The d s and y s are considered shift 

strains and account for the discontinuous through-thickness distributions of the strains due to the 
room-temperature bonding of the actuator to the laminate, and allow the strain in the actuator to be 
defined by reference-surface strains and curvatures. To explain, after the actuator is bonded to the 
laminate, the strains in the laminate would be due to cooling and due to bending from adding the 
actuator. The strains in the actuator would be due only to bending. As a result, the strain distributions 
through the thickness of the laminate and actuator would not be continuous. When the laminate- 
actuator combination is snapped to the second actuator-added shape, the laminate and actuator would 
both have extensional and bending strains; the strain distributions would continue to be discontinuous 
between the laminate and actuator. With actuation, it was assumed, in the spirit of the Kirchhoff 
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hypothesis, that the strain increments, in the laminate and actuator would be continuous through the 
thickness. As a result, the profdes of strains through the laminate and actuator would remain 
discontinuous. To account for this discontinuity, yet have one set of coefficients, i.e., c\ through c m, 
to define the strains within the laminate and the actuator due to the application of voltage, the shift of 
strain measures was necessary. The shift strains are thus given by 
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where the s'° s and y arc the initial strains and the k'° s are the initial curvatures, i.e., the strains 

and curvatures at the reference surface from Step II. These initial strains and curvatures are known. 
The total potential energy for the actuation portion of the model is thus 

n = n!+n 3 (19) 

The shapes, actuated or unactuated, of the laminate-actuator combination are determined by solving 
for ci through cio by equating to zero the first variation of the total potential energy of eq. 19. 
Stability is checked by examining at the second variation. 


COMPARISON OF MODEL WITH EXPERIMENTS 

For this study, a [0/90 ]t laminate, fabricated from AS4/3502 graphite-epoxy prepreg, was cured 
flat and then cooled to room temperature. A parameter study using the developed analytical model 
was conducted to determine the dimensions of the laminate, i.e., L x and L v , that could be snapped by 
the actuator. (The performance of the MFC actuator had been calibrated in a separate step using a 
cantilevered aluminum plate.) As determined by the parameter study, the laminate was cut to 150 by 
150 mm. The dimensions of the laminate and curvatures of the two cylindrical shapes were measured. 
An MFC actuator was centrally bonded to one side of the laminate to form the [0/90/0 MFC ] T laminate- 
actuator combination. 

The measurements of the laminate affected the inputs to the model as follows: When the 
curvatures of the laminate were measured, the major curvatures of the two stable shapes were not of 
the same magnitude. It was assumed that uneven resin bleed during manufacture was responsible for 
the unequal curvatures. In the model, it was found that by increasing the thickness of one layer and 
decreasing the thickness of the other layer a corresponding amount, the slightly different curvatures 
could be predicted quite accurately. Since the total thickness of the laminate was a measured quantity, 
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the assumption of the layers having an equal but opposite thickness change was justified. It was also 
found that, over time, the laminate lost some curvature, presumably due to moisture absoiption. The 
AT used in the model, in eq. 6, was reduced to account for this loss of curvature. The numerical 
values of the specific properties that were used in the analytical model are given in Table 1. 

To demonstrate the snap-through event, actuation experiments were performed. The experimental 
setup consisted of a voltage supply, voltage amplifier, strain gage amplifier/conditioner, voltmeter, a 
LabVIEW data acquisition system on a personal computer, and the laminate-actuator combination 
with a strain gage bonded to the side opposite the actuator. The output from the strain gage was used 
to determine when the snap through event occurred. To perform the experiment, the laminate-actuator 
combination was snapped, by hand, to the second actuator- added shape, fig. 2a. Since the analyses of 
the model assumed the edges of the laminate-actuator combination were free of any specified force 
resultants or displacements, the laminate-actuator combination was suspended, like a pendulum, by 
the strain gage wire when conducting the snap-through experiments. The voltage was increased from 
zero until the laminate-actuator combination snapped to a shape similar to that of fig. 2c; the voltage 
was then decreased to zero and the laminate-actuator combination assumed the shape given by fig. 2c. 
This procedure was repeated six times with good repeatability and an average snapping voltage of 
1695 V. 

With the adjustments to the layer thicknesses and AT, which were discussed above, the model 
was used to predict the behavior of the laminate as the voltage was varied. The values of C\, c 2 , c 3 , and 
c 4 were computed as a function of the voltage increase. These relationships are shown in fig. 3, the 
solid lines represent stable equilibrium shapes and the dashed lines represent the unstable equilibrium 
shapes. Recall, from eq. 2, that these coefficients describe the out-of-plane displacement of the 
laminate-actuator combination. As shown in fig. 3, for a given voltage level, each of these 
coefficients has either three solutions, or one solution. At voltage levels below 1262 V there are three 
solutions. For voltages greater than 1262 V, there is just one solution for each of the coefficients. 
Consider the laminate- actuator combination at zero applied voltage in the shape similar to fig. 2a, i.e., 
ci ~ -1 m 1 and c 3 ~ +95 m 2 , while c 2 and c 4 are both approximately zero. As the voltage is increased 
from zero, both ci and c 3 increase in magnitude, while c 2 and c 4 remain close to zero. At 1262 V, the 
solution exhibits limit-point behavior, and a further increase in voltage results in the solution jumping 
to a configuration with a large curvature in the v-direction and little curvature in the x-direction, i.e., 
Ci ~ 0, C 2 ~ -5.5 m ', c 3 ~ 0, c 4 ~ -45 m 2 . With a return of the voltage to zero, the room-temperature 
equilibrium shape similar to fig. 2c results. 

The two predicted room-temperature shapes of the laminate-actuator combination are shown 
along with the actual shapes in fig. 2. Though this figure illustrates only a qualitative comparison, 
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which is quite good, the quantitative comparison is good also. Considering the simplicity of the 
model, the predicted snapping voltage of 1262 V also compares fairly well with the snapping voltage 
measured in the experiments of 1 695 V. 


CONCLUDING REMARKS 

An interesting concept for moiphing structures has been presented. Modeling and experimental 
work was considered, and fairly good agreement between the model and experiment was found. 
Though this initial effort was successful, there are a number of directions in which this work can 
proceed. The most obvious direction would be to add another actuator to the other side of the 
laminate. This second actuator could be used allow the laminate to self-reset, i.e., snap from one 
shape to another and back to the original shape using only voltage applied to the actuators. The 
number of actuators on one or both sides, actuator dimensions and shapes, and actuator location can 
be varied to extend and build upon the particular problem studied. Changing the scale of the concept, 

i.e., changing laminate sidelength dimensions and number of layers, is also worth pursuing. 
Laminates with fiber angles other 0 and 90° could lead to shape changes that include twist curvature. 
In addition, more detailed analysis techniques, such as finite-element analysis, may be needed to 
capture some of more local effects that occur in this type of structure. Further details of the present 
effort are available in ref. 7. 
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Table 1 Material properties and other inputs to analytical model 


Property 

AS4/3502 
(Layer 1) 

AS4/3502 
(Layer 2) 

MFC Actuator 
(active portion) 

E\ (GPa) 

132.0 

128.0 

29 

E 2 (GPa) 

9.798 

9.608 

18 

Gn (GPa) 

5.112 

4.895 

5.3 

V\2 

0.2990 

0.3009 

0.28 

a, (1/°C) 

/ 1 /o r*\ 

-0.04 156x1 O' 6 

'll ^ 

5.538xl0' 9 

- 

Ct2 ( 1/ C) 

dn (//£/(kV/mm))* 

Z3. / /xlU 

Z4. /oXlU 

281 

dn (//A(kV/mm))* 

_ — — ^ 


-111 

thickness (m) 

133.8xl0' 6 

138.2xl0' 6 

290. xlO' 6 

L x (m) 

0.151 

8 

0.08573 

U (m) 

0.1518 

0.05715 

Axi (m) 



1.0668xl0' 3 

Ar(°c) 

-117 


* determined by calil 

bration with cantilevered aluminum plate 


FIGURE CAPTIONS 

Fig. 1 Geometry of laminate-actuator combination: (a) cross-section view, (b) top view 

Fig. 2 Experimental and predicted stable shapes of laminate-actuator combination with zero applied voltage 

Fig. 3 Coefficients C; through c 4 vs. applied voltage, showing existence of snap through 
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Fig. 3 Coefficients c; through c 4 vs. applied voltage, showing existence of snap through 
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